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Abstract. 

In evolutionary game dynamics, reproductive success increases with the 
performance in an evolutionary game. If strategy A performs better than strategy 
B, strategy A will spread in the population. Under stochastic dynamics, a single 
mutant will sooner or later take over the entire population or go extinct. We analyze 
the mean exit times (or average fixation times) associated with this process. We 
show analytically that these times depend on the payoff matrix of the game in an 
amazingly simple way under weak selection, i. e. strong stochasticity: The payoff 
difference Att is a linear function of the number of A individuals i, Att = ui + v. 
The unconditional mean exit time depends only on the constant term v. Given 
that a single A mutant takes over the population, the corresponding conditional 
mean exit time depends only on the density dependent term u. We demonstrate 
this finding for two commonly applied microscopic evolutionary processes. 
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1. Introduction 

Systems in which successful strategies spread by imitation or genetic reproduction 
can be described by evolutionary game theory. Such models are routinely analyzed 
in evolutionary biology, sociology, anthropology and economics. Recently, the 
apphcation of methods from statistical physics to these systems has lead to many 
important insights [H El |3], HI [5j. 

Traditionally, the dynamics is described by the rephcator equations, where the 
growth rate of a strategy is associated with its relative success compared with the 
population average PI7|. 

In the past years, research has focused on stochastic evolutionary game dynamics 
in finite populations [HliailDliniliaiiaiaiiaiiniinillllllSlEOlE^ In this 
context, a connection to the weak selection hmit of population genetics has been 
established [8j. Weak selection means that the payoff differences based on different 
strategic behavior in interactions represent only a small correction to otherwise 
random dynamics, similar to high temperature expansions in physics. Weak selection 
is considered as a relevant limit in biology, as most evolutionary changes are driven 
by small fitness differences [23]. Moreover, it allows analytical approximations that 
are often impossible when selective differences in payoffs are large P EH [25] . 

Most of the recent work that uses the weak selection approximation has been 
focusing on the probability that a certain strategy takes over. The time associated 
with this process has been calculated |26], but it received considerably less attention 
so far. Here, we present the weak selection corrections to the conditional and 
unconditional mean exit or fixation times in evolutionary 2x2 games with N players. 

The conditional average time to fixation is the expected time a single mutant 
needs to take over the population, given that such a takeover occurs at all. The 
unconditional average time of fixation ti is the expectation value for the time until 
the population is homogenous again after the arrival of a single mutant. This is 
regardless of wether the mutant type takes over the population or becomes extinct. 
Equivalently, the average fixation times for such one dimensional random walks can 
also be interpreted as mean first passage times or mean exit times [271 [281 [29] . 
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Throughout this paper, we use the payoff matrix 

A B 
A f a b 
Bye d 

An A player interacting with another A receives a. If it interacts with B, it obtains 
b. Similarly, B receives c from A and d from other i?'s. Thus, the average payoffs 
are 

, ■^ i — ^ N — i ^ 
-^W = ^« + ^& (2) 

i N — i — 1 

^^^^) = ]^"+ N-1 

A quantity that is of particular interest is the difference between the average payoffs, 

A-K{i) = TCAii) - TTBii) = ui + V, (4) 

where 

a + d - (b + c) 
iV-1 ' 
^ N{b-d)-{a-d) 

N -I ' ^ ' 

We show that under weak selection, the conditional time (t^) during which a single 
mutant takes over the whole population depends only on u (and, of course, on the 
population size). The unconditional time (ti) during which the mutant either takes 
over the population or reaches extinction depends only on v (and the population 
size). See Figure [T]for an illustration of the relevant quantities. 

Our manuscript is organized as follows: In Section [2j we introduce a particular 
evolutionary process for our analysis. Although our results are valid for a broader 
class of processes, we only present the full calculation for this evolutionary process. 
In Section [3| we recall the general form of fixation probabilities and times. We discuss 
neutral selection in Section [4] as a prerequisite to the weak selection expansion, which 
we explore in Section [5| In Section |6] we address the frequency dependent Moran 
process to underhne the generality of our findings. The consequences of our analytical 
results are discussed in Section [7l 
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Figure 1: Illustration of the most relevant quantities. We are interested in the evolutionary fate of 
a single A player. All quantities depend on the intensity of selection jS and the population size N. 
The payoff difference between A and B players is given by Att = ui + v, with i as the number of 
A players. Both the transition probabilities and T~ and the probability that a single A player 
takes over the population 0i depend on u and v. But for weak selection, /3 <C 1, the conditional 
time tf during which a single A player takes over a population of B players only depends on u, 
whereas the unconditional time ti until either A ot B has taken over the population only depends 
on V. 



2. Fermi process 

In a finite population of size N with two possible strategies A and B, the state of 
the system is characterized by the number of type A individuals i. In general, the 
dynamics is stochastic. In each time step, a randomly chosen individual evaluates 
its sucess. It compares this payoff with a second, randomly chosen individual. If 
this second individual has a higher payoff, the first one switches strategies with 
probability p > I- Otherwise, it switches with p < ^- We assume that the switching 
probability is given by the Fermi distribution. Its shape is controlled by the intensity 
of selection (3, which can be interpreted as an inverse temperature, 

r,± = I = ^ (7) 

In previous work |30l [311 IS^, there is a different strategy update procedure. The 
first individual switches to the second's strategy with probability pf. The second 
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individual can also switch to the first individual's strategy with probability 1 —pf- 
This yields a factor 2 in the transition probabilities (and, as we will become clear 
later, a factor ^ in the fixation times). This process also has a proper strong 
selection hmit, i. e. it is possible to examine (3 ^ oo. In this latter case we have 
pf — ^ 6(A7r(i)), where Q{x) is the step function. 

The population size is constant in time, in each time step the state of the system 
can at most change by one, i. e. from i to i — 1 or to 2 + 1. The transition probabilities 
Tf" to move from i to i ±1 are 

Tt = ^^pf- (8) 

The probability to stay in the current state is 1 — — . An important measure 
of where the system is more likely to move is their ratio, 

7. = ^=e-^^'^«. (9) 

i 

This is a quantity that describes the tendency to move from the state z to i =F 1, 
depending on whether 7^ ^ 1. Of course, > is required, which follows from 
/? < 00. The Tf^ and thus the 7^ are invariant under adding a value to each of the 
payoffs given in ([T]), whereas multiplying the payoff matrix with a factor A results in 
a change in the intensity of selection (3 = (3\. 

Let us now focus on weak selection, /5 -C 1. In this case we have 

pt^\± ^A7r(z). (10) 

Weak selection corresponds to high temperature in Fermi statistics. A Taylor 
expansion of the 7^ up to first order in (3 yields 7^ ^ 1 — /3A7r(i). In this case, 
the probability to move from i to i + 1 is very similar to the probability to move 
from i to i — 1. Weak selection links the Fermi process to a variety of birth death 
processes, cf. P 133]. 



3. Fixation probabilities and fixation times 

From equation ^ it follows that the two pure states all A or all B are absorbing, 
Tq = = 0. In a finite population, we can calculate the probability 0j that the 
system will fixate to the pure state all A, starting with the mixed state i. Obviously, 
we have 0o = and (pN = 1. For < i < A^, there is a balance equation for the 
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fixation probabilities, (pi = T- + (1 — — T. )0j + T-+0j+i. This recursion leads 
to an expression for the fixation probabilities in terms of the 7^ |3ll |35l 136] . 

, _ i + Ei=\nti7i 



1 I ^N-l T-{k ■ 

1 + Ek=i Ui=i ii 
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which is valid for any birth death process. 

For the Fermi process, the exact equation ^ simplifies matters in an elegant 



way because the products in equation (11) can be solved, 

k 



= exp \ -(3 5^ A7r(/) = exp {-/? [fc^^ + k{^+v)\]. (12) 



1=1 V 1=1 



Hence, equation (11) simplifies to 



(13) 



l + Ef=Yexp{-/5[fc2| + A;(| + ^;)]}- 

For large N , the sums in equation ( [T3| can be approximated by integrals, which 
yields a closed expression for the probabilities 0j [33, EI]- 

General expressions for the unconditional and conditional mean exit times 
or average times of fixation, ti and t^, are well known, especially for simple, 
translational invariant random walks [211 EH ES]- A complete derivation for the 
average times of fixation in finite systems without translational invariance can be 
found in [ISlESlEn]. 

In the following, we will focus on the fixation of a single A mutant in a population 
of B. Accordingly, the unconditional and conditional fixation times read 

A'— 1 k k 
k=l 1=1 ' m=l+l 

and 

N-l k , k 

'i'^EElt n (15) 

k=l 1=1 I m=l+l 

respectively. Time is measured in elementary time steps here. Thus, in each time 
step one reproductive event occurs. In biological contexts, it is often more convenient 
to measure time in generations, such that each individual reproduces once per 
generation on average. Time in generations is obtained by dividing the number 
of time steps by the population size N. It is well known that the variance of the 
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exit times under weak selection can be large [32], which has important biomedical 
implications [lO]. Nonetheless, here we concentrate on the expectation values and 
do not address the distribution of the exit times. 



4. Neutral selection 



An important reference case is neutral selection, which results from vanishing 
selection intensity (3 = [H]. Neutral selection is a very general limit, which is 
typically not affected by the details of the evolutionary process. For neutral selection 
we have 7^ = 1, that is = in any state i. However, we still have 7^ 
for i 7^ j, although the system is symmetric, T^^ = T^_j. This is a difference to the 
simple random walk in one dimension, which is invariant with respect to translation 

m- 

For the Fermi process, the neutral transition probabilities are 

1 i N-i 



Tf 



1. From equation (11), it is thus clear that 



^=0 2N N 
We have = T~, which leads to 7, 
the probability of fixation to A is given by the initial abundance of A, 

i 

For the neutral unconditional time of fixation ti we get 

1 'v^^v^ 2N^ 



(16) 



13=0 



N-l k 

N 



2NH 



N-l- 



(17) 



(18) 



k=l 1=1 ^ ' 

Details for this calculation can be found in [Appendix A[ We introduced the shorthand 
notation for the harmonic numbers H^^x = YliJi^ which diverge logarithmically 
with N. In the same way we can solve 

For neutral selection, the conditional average time of fixation of a single mutant 
diverges quadratically with the system size. 
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5. Weak Selection 

In this section we will calculate the linear corrections of the mean exit times or 
fixation times ti, and under weak selection, /? -C 1. Of course, all weak selection 
approximations are valid only if the term linear in f3 is small compared to the constant 
term. 

The fixation probabilities for small f3 are 
which has been derived for a variety of evolutionary processes before [HI HOl [El EHl 



Next, we address the weak selection approximation of the fixation times. The 
expectation value of the unconditional fixation time of a single A mutant in a 



population of -B is in general given by the exact equation (14). With the the transition 
and fixation probabilities of the Fermi process, the unconditional fixation time of 
absorption at any boundary simplifies to 



k= 



X exp i -(3 ^ A7r(m) I . (21) 



m=l+l 



The weak selection approximation takes the remarkably simple form (see Appendix 
B] for details) 

ti^2NHN-i+vN{N -I-Hn-i) (3, (22) 

with V given in ([g]). Thus, ti depends only on the constant term of the payoff 
difference. For large iV, this yields v ^ b — d. That is, for large populations under 
weak selection the linear correction of the average fixation time only depends on 
the advantage (or disadvantage) of the A mutants in the resident population. For 
b > d, invasion of A mutants is likely and slows down the time until the population 
is homogeneous again. For > 6, it is difficult for A to invade a B population and 
extinction of the mutants is faster than in the neutral case. Note that the payoff 
entries a and c have no infiuence on the unconditional fixation time under weak 
selection corrections. Since fixation is unlikely for weak selection (the probability of 
fixation of a single A mutant is approximately N~^), the unconditional fixation time 
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is dominated by the fixation to B. In this case, it is enough to discuss the invasion 
of A mutants. 

Next, we address the average time to fixation given that the A mutant takes 
over the population. With the general result (15) the Fermi processes conditional 
fixation time to all A reads 

k=l 1=1 ^ ' 

xexp|-/3 ^ A7r(m)l. (23) 

L m=l+l ) 

Its linear approximation turns out to be dependent on the payoff's in a very simple 
way as well, 

jV^ + - 6 

^2N{N -l)-uN{N -I) -(3, (24) 

18 



with u{N — l) = a — h — c + d. The detailed calculation can be found in Appendix 



[B| Since during the fixation process all payoffs are of importance, it is obvious that 
they all enter here. For example, when it is easy to invade because few mutants 
have an advantage {h > d), but difficult to reach fixation because mutants are 
disadvantageous once they are frequent (c > a), we have u < and the conditional 
time to fixation is larger than neutral. In the last section, we discuss special classes of 
games to show that, under weak selection, the conditional mean exit times of fixation 
(or absorption) do not always follow the intuition based on the payoff matrix ([T|. 



6. Frequency dependent Moran process 

In this section we address the generality of the previous findings discussing an 
alternative evolutionary process. The first model that connects payoffs from a 2 x 2 
game to reproductive fitness using a weak selection approach in finite populations 
is the frequency dependent Moran process [H |9]. In this process, an individual is 
chosen for reproduction with probability proportional to its fitness f{i). The offspring 
replaces a randomly chosen individual. The average payoffs ([2]) and ^ are mapped 
to the fitness such that /yi(«) = I — /3 + (3 7r^(i) and /b(«) = 1 — /3 + (3 vrB(«), where 
the selcetion intensity /3 > is so small that /a(«) > and fsii) > 0. The transition 
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probabilities of the standard Moran process read 
+ _ N - I 

' rf^{r) + [N-i)fs{^ AT ' ^ ^ 

^ {N-i)fB{i) i_ . . 

' ^f^{{) + {N~^)M^)N■ ^ > 

Although these transition probabilities are different from those of the Fermi process, 
they also yield 7^ ^ 1 - A7r(z) and nm=/+i 7m ~ 1-/3 Y.t,=i+i A7r(m) for 
weak selection, /5 <C 1. Thus, the weak selection approximations of the fixation 
probabilities (pi of the Moran process and the Fermi process are identical, see equation 



(20). But the weak selection approximation of the transition probabilities are not 
identical, which leads consequently to different mean exit times. Nevertheless, the 
results have the same, remarkably simple connection to the payoff matrix ([T]). The 
mean exit times or fixation times of the frequency dependent Moran process are 

ti ^ NH^^, + ^ _ (iv + 1 - 2Hn) /?, (27) 

tt^NiN-D-u'^^i^'^^^p. (28) 

Qualitatively, the dependence on the payoff matrix via u and v is the same as for the 
Fermi process. Their calculation is analogous to the findings of the previous section, 
details can be found in [Appendix B[ Note that, comparing with the Fermi process. 



there is a factor of 2 missing in the neutral terms. However, this can be avoided by 
rescahng the transition probabilities, without changing the properties of the different 
processes. 



7. Discussion 



Finally, let us discuss the imphcations of our results for general 2x2 games. While 
we concentrate on the Fermi process here, the discussion is equally vahd for the 
frequency dependent Moran process. An important question is whether the hnear 
correction for weak selection is compatible with the general features of the game and 
the known asymptotic behavior for large N of the mean exit or fixation times derived 
by Antal and Scheuring [261 • Clearly, this depends on the payoff matrix of the 2x2 
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game, 



(29) 

as the payoffs enter tfie first exit times of absorption linearly. To analyze the 
difference to the neutral case we consider the rescaled average times of fixation, 
Ti{/3) = ti{/3)/ti{0) and t^{/3) = tf{/3)/tf{0). The rescaled unconditional fixation 
time reads 

lN(b-d)-a + d fN-1 \ ^ 

Accordingly, the rescaled conditional fixation time for absorption at all A is 
, ^ a-b-c + d + 
^^-^ N-1 36 ^^^^ 

Note that for population sizes N > 2 and sufficiently small /3, we always have ti{0) < 
ti{0). In other words, the average time until the A individual has reached fixation 
or gone extinct is smaller than the conditional average time until the A individual 
has reached fixation. For P ^ oc, the process follows deterministically the intensity 
of selection and thus both fixation times may coincide, ti{f3 ^ oo) ^ ti{(3 ^ oo). 
This ordering of the fixation times is blurred by our rescaling, as we focus only on 
the change relative to the neutral case. 

In the following, we discuss these two expressions for the three generic types of 
2x2 games, namely dominance oi A {a > c and h > d), coexistence of A and B 
{a < b and c> d) and a coordination game (a > c and b < d). 



7.1. Dominance of A. 

Consider a game where strategy A is always dominant, i. e. it obtains a larger payoff 
than B, regardless of the fraction of A in the population. This is the case for a > c 
and b > d. One special case is the Prisoner's Dilemma with b > d > a > c. 
The interesting feature of this game is that the social optimum d is not the Xash 
equilibrium, which is a. For neutral selection, a single A individual goes extinct with 
probability 1 — N~^. Thus, the unconditional fixation time ri is dominated by the 
extinction of A. Since strategy A is favored by selection, increasing the intensity 
of selection decreases the probability of the extinction of A. Since fixation takes 



Fixation times in evolutionary games under weak selection 



12 



at least N — 1 time steps, ti increases with increasing intensity of selection (3. For 
large N, this is obvious from our equation (30), because in this case the quantity 
N{h—d) — a+d is positive. However, once extinction of A becomes unhkely, increasing 
/? further will lead to a decrease of ri. 

The discussion of the conditional fixation time is not as straightforward, 
because the sign oia — h — c + d can be positive or negative. The sign of this 
quantity is also decisive for the evolutionary dynamics in other contexts, see e.g. 
[33]. When the advantage of an A individual is initially large and decreases with 
the abundance of A (a — c > 6 — > 0), then the sign of a — h — c + d is positive 
and Ti decreases with increasing intensity of selection. But when the advantage 
of strategy A decreases with the number of A individuals (6 — (i>a — c>0), 
then increases with increasing intensity of selection. However, this apparently 
counterintuitive phenomenon (after all, A dominates B) can only be observed for 
weak selection. For strong selection, decreases again. These results are compatible 
with the observation that the conditional fixation time scales as A^lnA^ for large 
[26] . In Figure [2] (a) we show a numerical example for the rescaled average times. We 
include averages from numerical simulations of the evolutionary process, our linear 
approximation as well as the exact result that can be obtained from dividing equation 
(14) by (18) and equation (15) by (19), respectively. The payoff matrix is chosen 
such that a + d > h + c, which means that with increasing intensity of selection 
decreases and ri increases. 



7.2. Coexistence of A and B. 

As a second class, we consider games in which B is the best reply io A {c> a), but 
A is the best reply to B {b > d). Important examples for such games are the Hawk- 
Dove game [35 the Snowdrift game [3S]- For infinite populations, the rephcator 
dynamics predicts a stable coexistence oi A and B. In finite populations, the system 
typically fluctuates around that point until eventually, fluctuations lead to absorption 
in one the boundaries [301 131] • Consequently, the conditional fixation times increase 
exponentially with the population size [26]. Since a — b — c + d is negative, we 
also have an increase of with the selection intensity for weak selection. Further, 
N{b ~ d) — a + d is positive in large populations, such that also ti increases with 
the selection intensity. Figure [2] (b) shows that the divergence of the exact results is 
faster than the linear approximation even for weak selection. 



Fixation times in evolutionary games under weak selection 



13 



7.3. Coordination games. 

Finally, let us discuss coordination games in which a > c and b < d. In these games, 
A is the best reply to A and B is the best reply to B. The replicator equation of 
such systems exhibits a bistability: If the fraction oi A individuals is sufficiently high 
in the beginning, the A individuals will reach fixation. Otherwise, B individuals 
will take over the system. The stronger the intensity of selection, the less likely it 
is that a single A individual can take over a B population. Consequently, ri should 
decrease with /3. This also follows from our weak selection approximation: In large 
populations, N{b — d) — a + d is negative and thus ri decreases with the intensity of 



selection, see equation (30). Perhaps less intuitive, also rf decreases with /?, which 



results from a — b — c + d> 0, cf. (30). However, this is again consistent with the 
observation that scales as A^lniV in large populations. Although the fixation 
probability of a single A decreases with /?, if such an event occurs, it is faster than 
in the neutral case. A numerical example for this behavior is shown in Figure [2] (c). 

The numerical examples indicate that the convergence radius of our weak 
selection expansion is of the order of A^~^, which is also known for many systems in 
population genetics. Although might appear small, this kind of weak selection 
is the most relevant limit in evolutionary biology, as evolutionary change is typically 
only connected with small selective differences. We stress that we have made no 
assumptions on the population size, such that our results are vahd for arbitrary A^. 

Our approach shows under which circumstances the general features of the game 
are refiected in the fixation times under weak selection. Although the weak selection 
expansion of the mean exit or fixation times is technically rather tedious, the resulting 
asymptotic behavior shows remarkable simplicity. 
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Figure 2: Expectation values of the rescaled fixation times starting with a single A mutant 
in a population of i? as a function of the selection intensity f3. Full lines show the normalized 



exact solution originating from the exact results (21 1 and (23 1. Colored dashed lines are the linear 



approximations (30 1 and (31 1. Symbols show the results from simulations based on 10'' realizations, 



which agree nicely with the exact results. Diamonds are for the unconditional averages, circles are 
for the conditional averages. On the right hand side, we show the payoff matrices of the three games 
and illustrate the direction of selection in these games, (a) In a game with dominance of strategy 
A, the unconditional fixation time increases with the intensity of selection, but the conditional 
fixation time decreases, (b) For games with stable coexistence, both fixation times increase with 
the intensity of selection, (c) For coordination games, the two fixation times become shorter when 
the intensity of selection is increased. In all examples, the population size is = 100. 
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Appendix A. Finite double sums 

Here, we collect some helpful calculations for double sums as they appear in the 
mean exit times. An important observation is 

N~l k J. i-1 „ N-1 

k=i 1=1 1=1 l=i 

for any function // < oo and Z = 1, . . . , — 1. This can be seen by writing the left 
hand side term by term, i. e. 

N-l k „ „ J. 

\ ^ \ ^ Jl _ Jl Ji 

N -I ~ N -\ N -i 

k=i l=\ 

_|_ /l _|_ _|_ _j_ 



A^-1 N -i N 

+ ... 

' +••• + #- + •••+ (A.2) 



A^-1 N -i A^-(A^-l) 



N-l 

1=1 

i-1 J. N-l 

-('v-)E]^ + E/.- 

1=1 1=1 

For the case i = 1 the result is especially simple, since the first sum of the right hand 
side of equation (A.l) vanishes. This case is of special interest for the computation 
of under neutral selection with fi = l and for ti with // = l/l. 

Another finding for double sums with M eN and two bounded functions fk and 

9i is 

M k MM 
k=l 1=1 1=1 k=l 

This becomes clear by resorting the terms again, 

M k 

^^fk9l = flQl + f2{,gi +92) + ■■■ + fMigi + 92 + ■■■+ qm) 
k=l 1=1 
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M M 



Appendix B. Fixation times under weak selection 

Here, we calculate the linear corrections of the mean exit times ti and for the 



Fermi process in detail, compare equations (21) and (23). We aim at finding these 
times for weak selection, e.g. 

' ■ ^ (B.l) 



dp 



/3=0 



The first term follows directly from the calculation in Appendix A\ see equation (|18D 

d 4 



Our goal here is to compute the linear term 

N~l k r 



9/3 



d_ 
dp 



EE 

k=l 1=1 

N-1 k 

EE 



+ 



/3=0 

d 1 



T+ dp ' ^'dpT+ 



f3=0 



k=l 1=1 



(B.2) 



I E A' 

' m=l+l 
d 1^ 11™=/+! 7n 

For the fixation probability under weak selection and with A7r(/) = ul + v, we have 



where we applied 11™=;+! 7m | ^ o ^ ""^ 



(3=0 
13=0 



J2m=l+1 ^ 



■K[m). 



dp 



(3=0 



J_ (A + Ow + St; 

A^ ' 6 



(B.3) 



The weak selection approximation of the inverse of the transition probability T^, 
compare equation (Isj), yields 



A2 



d 1 



dPT+\., l{N-l) 



(ul + v). 



(B.4) 



Thisleads to 



d_ 
dp 



(3=0 



N-1 k 

EE 

k=l 1=1 

N~l k 



(A- 1)((A + l)u + 3w) 2A2 



6A 



1{N - I) 



EE A, 



^^A/(A-/) 



{ul + v) 
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=1 1=1 ^ ' m=l+/ 



N-\ k ^ 

-EE AT 



(B.5) 



While the first two double sums can be solved with the help of [Appendix A[ the third 



term is more complicated. For this more tedious calculation, we refer to [Appendix 
[Cj Eventually, the solution of the double and triple sums leads to 



d_ 
dp 



- A^(A^ - l)u - NHn^iv 
vN{N- 



Hn-1 



u 



(B.6) 



where the last step is elementary. Combining this with equation (18) leads finally to 
the unconditional mean exit time under weak selection, equation (22). 
For the conditional fixation time t^, the linear term 



a_ 4.A 

dl3 ''1 



reads 



J/3=0 



_d_ 

d(3 



N-l k 



13=0 



EE 

k=l 1=1 
N-l k 

EE 

k=l 1=1 



1 90; 



+ 



d 1 



dp ^ dpT+ 



m=l+l 



/3=0 



13=0 



(B.7) 



The linear term of the weak selection 



The only difference compared to the unconditional fixation time, equation ([B.2[), 
is the fixation probability (pi instead of (pi. 
expansion of 0; is given in equation (B.3). This yields 

N~l k 



d_ 
'dp 



ti 



(3=0 



EE 

k=l 1=1 

N~l k 



1{N - 1){{N + l)u + 3v) 2N^ 



EE 

k=l 1= 
N-l k 

EE 

k=l 1-- 



QN l{N-l) 

[ul + v) 

k 

{uni + v). 



I AT^ 
^ Nl{N -I) 

^ Nl{N-l) 



(B.8) 



m=l+l 



Again, the first two double sums can be solved using the results from [Appendix 



\K\ The third term follows from a calculation which is similar to Appendix C , but 
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simpler. This last term reduces to 



Af-l k 

EE 



(A^-2)(A^- 1) 



(um + v) 



m=l+l 



N 



{{5N + 3)u + 9v) . 



(B.9) 



Finally, combining the three terms again results in 



N^{N -1) 



{{AN + l)u + 9v) 



13=0 



N{N-l) 



-{N 



{Nu + 2v) 
N{N-1 



{{5N + 3)u + 9v) 



uN{N-l) 



(B.IO) 



In combination with equation (19), this results in the conditional mean exit time 
under weak selection, equation (24). 

For completeness, we briefly repeat this calculation for the mean exit times of 
the frequency dependent Moran process. With the transition probabilities (25) and 
(26), the fixation probabilities under weak selection are identical to those of the Fermi 
process, see equation (20). However, the inverse transition probabflity is different in 
the weak selection regime, i. e. the linear correction is 

'd_J_ 

Hence, for the unconditional mean exit time we have the same starting equation 
(B.2). But with equation (|B.11|) this gives 

d 



-A.(0 = -N— 



fB.lll 



d(3 



(5=0 



N-1 k 

EE 

k=l 1=1 



(A^- l)((iV + l)u + 3^;) A^^ 



N-1 k 

EE 

k=l 1=1 

N-1 k 

k 



N 



N 



6N 
ul + V 



1{N - I) 



I 



(B.12) 



m=l+l 
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which differs from equation (B.5) only in the second double sum. With the previous 



findings for the Fermi processes times the required calculation is straightforward and 
results in 



d_ 

d(3 



13=0 



N 

v-{N + 1-2Hn) 



(B.13) 



That is, this hnear correction has a different dependence on the system size N. 
For the conditional mean exit time the situation is similar. In difference to equation 



(B.8), the linear correction reads 

N~l k 



d_ 
dp 



1{N -1){{N + l)u + 3v) N'^ 



J/3=0 



EE 

k=l 1=1 

N-l k , , , 

^ ^ N 

k=l 1=1 

N-l k 



6N 



I 



1{N -I) 



; r (um + v). 

^ ^ NUN -I) ^ ^ ' 

k=l 1=1 ^ ' m=l+l 



This the leads to 



d_ 



ti 



13=0 



Ar2 

-u—{N'-3N + 2) 



(B.14) 



fB.15) 



for the linear correction of the conditional mean exit times of the frequency dependent 
Moran process. 

Appendix C. Finite triple sum 



Here, we calculate the triple sum from Appendix B that require some additional 
steps. Our goal is to solve 

N—l k k 

- = EE^E^'M' (C.I) 

fc=l 1=1 ^ ' m=l+« 

For the sum over payoff differences, we have 

k k 

A7r(m)= Y,[um^v) = ]^-ju (C.2) 

m=l+i m=l+l 

where we introduced the function 



u 



fm = rn{m + 1) — + m w , 



(C.3) 
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N 



1 I 1 
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EE 



fk 



N-l k 

EE 



f, 



N-l k 



N-l k 



fk 




EE EE 



' l{N-l) I ' N-P 



+ 



(C.4) 



We solve each part separately, starting with the last one. For K4, we obtain with 
equation (A.l) from [Appendix A 

N-l 

(C.5) 



N-l k „ 

'^■-EE^-E/' 

k=l 1=1 k=l 



N-l 

k = -^—N{{N+l)u + 3v). 



6 



The second last term, i^a, is a sum over a linear function and can be treated with 
any table of elementary sums, e. g. [38] , 



7V-1 k 



^<-3=EEt = EE(('+i 



N-l k 



u 



k=l 1=1 

N-l 



k=l 1=1 



12 



N{{N + A)u + Qv). 



(C.6) 



The remaining two terms require more effort. Both terms, Ki and K2 have the same 

of 

fk 

k=l 1=1 

and 



structure regarding functions of k and I. Using equation (A. 3), we have 

7V-1 k 

'^■^"EE 



N-l 



N-l ^ Af-1 



N-l 

1=1 k=l 



(C.7) 



7V-1 k 



k=l 1=1 1=1 k=l 

Hence, we first have to compute the sum X^fc^^ /fc? which reduces to the solution of 
elementary sums. 



Af-1 Af-1 

EjEA- 



k=l 



k=l 



Af-1 
k=l 



U 



E^ 

k=l 



Fixation times in evolutionary games under weak selection 

^ {N'' + Nl + l^ - l)u + 3(iV + / - l)v) 



21 



6 

^fN' + m + l'~l)u + ^{N + l-l)v. (C.9) 



Thus, solving equations (C.7) and (C.8) simplifies to solving the elementary sums 
J2i=i^ /'^ with s = 0, 1, 2, compare [48]. With this, we have 



AT 1 AT , Q 

K, = ^—J2iiN + l)u + 3v) + — 5^ Z + ttE^ 



6 

A^- 1 



1=1 



1=1 



1=1 



For Ki, we obtain 
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(C.IO) 



K 



N{N - 1)(A^ + l)u + 3N{N - l)v 
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N~l 
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((iV + l)„ + 3i;) 



6 

N-l 
36 



{N -2){{2N + 3)u + 9v) 



Summing up the terms, a = {Ki + K2 — K3 — K4)/N, finally yields the result 

N-l 

{{{N +1)u + 3v)Hn-i-3{u + v)). (C.ll) 



6 



Again, if„ = X]r=i ^1^ harmonic numbers. In equation (B.9), the reasoning 

is very similar, but only terms of the structure of K2 and appear. 
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